Current challenges and future of agricultural genomes to phenomes in the USA.

Dramatic improvements in measuring genetic variation across agriculturally relevant populations (genomics) must be matched by improvements in identifying and measuring relevant trait variation in such populations across many environments (phenomics). Identifying the most critical opportunities and challenges in genome to phenome (G2P) research is the focus of this paper. Previously (Genome Biol, 23(1):1-11, 2022), we laid out how Agricultural Genome to Phenome Initiative (AG2PI) will coordinate activities with USA federal government agencies expand public-private partnerships, and engage with external stakeholders to achieve a shared vision of future the AG2PI. Acting on this latter step, AG2PI organized the "Thinking Big: Visualizing the Future of AG2PI" two-day workshop held September 9-10, 2022, in Ames, Iowa, co-hosted with the United State Department of Agriculture's National Institute of Food and Agriculture (USDA NIFA). During the meeting, attendees were asked to use their experience and curiosity to review the current status of agricultural genome to phenome (AG2P) work and envision the future of the AG2P field. The topic summaries composing this paper are distilled from two 1.5-h small group discussions. Challenges and solutions identified across multiple topics at the workshop were explored. We end our discussion with a vision for the future of agricultural progress, identifying two areas of innovation needed: (1) innovate in genetic improvement methods development and evaluation and (2) innovate in agricultural research processes to solve societal problems. To address these needs, we then provide six specific goals that we recommend be implemented immediately in support of advancing AG2P research.

Integrative profiling of gene expression and chromatin accessibility elucidates specific transcriptional networks in porcine neutrophils.

Neutrophils are vital components of the immune system for limiting the invasion and proliferation of pathogens in the body. Surprisingly, the functional annotation of porcine neutrophils is still limited. The transcriptomic and epigenetic assessment of porcine neutrophils from healthy pigs was performed by bulk RNA sequencing and transposase accessible chromatin sequencing (ATAC-seq). First, we sequenced and compared the transcriptome of porcine neutrophils with eight other immune cell transcriptomes to identify a neutrophil-enriched gene list within a detected neutrophil co-expression module. Second, we used ATAC-seq analysis to report for the first time the genome-wide chromatin accessible regions of porcine neutrophils. A combined analysis using both transcriptomic and chromatin accessibility data further defined the neutrophil co-expression network controlled by transcription factors likely important for neutrophil lineage commitment and function. We identified chromatin accessible regions around promoters of neutrophil-specific genes that were predicted to be bound by neutrophil-specific transcription factors. Additionally, published DNA methylation data from porcine immune cells including neutrophils were used to link low DNA methylation patterns to accessible chromatin regions and genes with highly enriched expression in porcine neutrophils. In summary, our data provides the first integrative analysis of the accessible chromatin regions and transcriptional status of porcine neutrophils, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the utility of chromatin accessible regions to identify and enrich our understanding of transcriptional networks in a cell type such as neutrophils.

Toward a data infrastructure for the Plant Cell Atlas.

We review how a data infrastructure for the Plant Cell Atlas might be built using existing infrastructure and platforms. The Human Cell Atlas has developed an extensive infrastructure for human and mouse single cell data, while the European Bioinformatics Institute has developed a Single Cell Expression Atlas, that currently houses several plant data sets. We discuss issues related to appropriate ontologies for describing a plant single cell experiment. We imagine how such an infrastructure will enable biologists and data scientists to glean new insights into plant biology in the coming decades, as long as such data are made accessible to the community in an open manner.

Assessment of DNA methylation in porcine immune cells reveals novel regulatory elements associated with cell-specific gene expression and immune capacity traits.

BACKGROUND: Genetics studies in the porcine immune system have enhanced selection practices for disease resistance phenotypes and increased the efficacy of porcine models in biomedical research; however limited functional annotation of the porcine immunome has hindered progress on both fronts. Among epigenetic mechanisms that regulate gene expression, DNA methylation is the most ubiquitous modification made to the DNA molecule and influences transcription factor binding as well as gene and phenotype expression. Human and mouse DNA methylation studies have improved mapping of regulatory elements in these species, but comparable studies in the pig have been limited in scope. RESULTS: We performed whole-genome bisulfite sequencing to assess DNA methylation patterns in nine pig immune cell populations: CD21+ and CD21- B cells, four T cell fractions (CD4+, CD8+, CD8+CD4+, and SWC6γδ+), natural killer and myeloid cells, and neutrophils. We identified 54,391 cell differentially methylated regions (cDMRs), and clustering by cDMR methylation rate grouped samples by cell lineage. 32,737 cDMRs were classified as cell lowly methylated regions (cLMRs) in at least one cell type, and cLMRs were broadly enriched in genes and regions of intermediate CpG density. We observed strong correlations between differential methylation and expression across immune cell populations, with cell-specific low methylation disproportionately impacting genes exhibiting enriched gene expression in the same cell type. Motif analysis of cLMRs revealed cell type-specific enrichment of transcription factor binding motifs, indicating that cell-specific methylation patterns may influence accessibility by trans-acting factors. Lastly, cDMRs were enriched for immune capacity GWAS SNPs, and many such overlaps occurred within genes known to influence immune cell development and function (CD8B, NDRG1). CONCLUSION: Our DNA methylation data improve functional annotation of the porcine genome through characterization of epigenomic regulatory patterns that contribute to immune cell identity and function, and increase the potential for identifying mechanistic links between genotype and phenotype.

Intestinal single-cell atlas reveals novel lymphocytes in pigs with similarities to human cells.

Lymphocytes can heavily influence intestinal health, but resolving intestinal lymphocyte function is challenging as the intestine contains a vastly heterogeneous mixture of cells. Pigs are an advantageous biomedical model, but deeper understanding of intestinal lymphocytes is warranted to improve model utility. Twenty-six cell types were identified in the porcine ileum by single-cell RNA sequencing and further compared with cells in human and murine ileum. Though general consensus of cell subsets across species was revealed, some porcine-specific lymphocyte subsets were identified. Differential tissue dissection and in situ analyses conferred spatial context, revealing similar locations of lymphocyte subsets in Peyer's patches and epithelium in pig-to-human comparisons. Like humans, activated and effector lymphocytes were abundant in the ileum but not periphery of pigs, suggesting tissue-specific and/or activation-associated gene expression. Gene signatures for peripheral and ileal innate lymphoid cells newly discovered in pigs were defined and highlighted similarities to human innate lymphoid cells. Overall, we reveal novel lymphocyte subsets in pigs and highlight utility of pigs for intestinal research applications.

Swine models for translational oncological research: an evolving landscape and regulatory considerations.

Swine biomedical models have been gaining in popularity over the last decade, particularly for applications in oncology research. Swine models for cancer research include pigs that have severe combined immunodeficiency for xenotransplantation studies, genetically modified swine models which are capable of developing tumors in vivo, as well as normal immunocompetent pigs. In recent years, there has been a low success rate for the approval of new oncological therapeutics in clinical trials. The two leading reasons for these failures are either due to toxicity and safety issues or lack of efficacy. As all therapeutics must be tested within animal models prior to clinical testing, there are opportunities to expand the ability to assess efficacy and toxicity profiles within the preclinical testing phases of new therapeutics. Most preclinical in vivo testing is performed in mice, canines, and non-human primates. However, swine models are an alternative large animal model for cancer research with similarity to human size, genetics, and physiology. Additionally, tumorigenesis pathways are similar between human and pigs in that similar driver mutations are required for transformation. Due to their larger size, the development of orthotopic tumors is easier than in smaller rodent models; additionally, porcine models can be harnessed for testing of new interventional devices and radiological/surgical approaches as well. Taken together, swine are a feasible option for preclinical therapeutic and device testing. The goals of this resource are to provide a broad overview on regulatory processes required for new therapeutics and devices for use in the clinic, cross-species differences in oncological therapeutic responses, as well as to provide an overview of swine oncology models that have been developed that could be used for preclinical testing to fulfill regulatory requirements.

Pig genome functional annotation enhances the biological interpretation of complex traits and human disease.

The functional annotation of livestock genomes is crucial for understanding the molecular mechanisms that underpin complex traits of economic importance, adaptive evolution and comparative genomics. Here, we provide the most comprehensive catalogue to date of regulatory elements in the pig (Sus scrofa) by integrating 223 epigenomic and transcriptomic data sets, representing 14 biologically important tissues. We systematically describe the dynamic epigenetic landscape across tissues by functionally annotating 15 different chromatin states and defining their tissue-specific regulatory activities. We demonstrate that genomic variants associated with complex traits and adaptive evolution in pig are significantly enriched in active promoters and enhancers. Furthermore, we reveal distinct tissue-specific regulatory selection between Asian and European pig domestication processes. Compared with human and mouse epigenomes, we show that porcine regulatory elements are more conserved in DNA sequence, under both rapid and slow evolution, than those under neutral evolution across pig, mouse, and human. Finally, we provide biological insights on tissue-specific regulatory conservation, and by integrating 47 human genome-wide association studies, we demonstrate that, depending on the traits, mouse or pig might be more appropriate biomedical models for different complex traits and diseases.

Associations of natural variation in the CD163 and other candidate genes on host response of nursery pigs to porcine reproductive and respiratory syndrome virus infection.

Pigs with complete resistance to porcine reproductive and respiratory syndrome (PRRS) virus (PRRSV) have been produced by genetically knocking out the CD163 gene that encodes a receptor of the PRRSV for entry into macrophages. The objectives of this study were to evaluate associations of naturally occurring single nucleotide polymorphisms (SNPs) in the CD163 gene and in three other candidate genes (CD169, RGS16, and TRAF1) with host response to PRRSV-only infection and to PRRS vaccination and PRRSV/porcine circovirus 2b (PCV2b) coinfection. SNPs in the CD163 gene were not included on SNP genotyping panels that were used for previous genome-wide association analyses of these data. An additional objective was to identify the potential genetic interaction of variants at these four candidate genes with a mutation in the GBP5 gene that was previously identified to be associated with host response to PRRSV infection. Finally, the association of SNPs with expression level of the nearby gene was tested. Several SNPs in the CD163, CD169, and RGS16 genes were significantly associated with host response under PRRSV-only and/or PRRSV/PCV2b coinfection. The effects of all SNPs that were significant in the PRRSV-only infection trials depend on genetic background. The effects of some SNPs in the CD163, CD169, and RGS16 genes depend on genotype at the putative causative mutation in the GBP5 gene, which indicates a potential biological interaction of these genes with GBP5. In addition, genome-wide association results for the PRRSV-only infection trials revealed that SNPs located in the CDK5RAP2 or MEGF9 genes, near the TRAF1 gene, had suggestive effects on PRRS viral load, which indicates that these SNPs might contribute to PRRSV neuropathogenesis. In conclusion, natural genetic variants in the CD163, CD169, and RGS16 genes are associated with resistance to PRRSV and/or PCV2b infection and appear to interact with the resistance quantitative trait locus in the GBP5 gene. The identified SNPs can be used to select for increased natural resistance to PRRSV and/or PRRSV-PCV2b coinfection.

Quantitative analysis of the blood transcriptome of young healthy pigs and its relationship with subsequent disease resilience.

BACKGROUND: Disease resilience, which is the ability of an animal to maintain performance under disease, is important for pigs in commercial herds, where they are exposed to various pathogens. Our objective was to investigate population-level gene expression profiles in the blood of 912 healthy F1 barrows at ~ 27 days of age for associations with performance and health before and after their exposure to a natural polymicrobial disease challenge at ~ 43 days of age. RESULTS: Most significant (q < 0.20) associations of the level of expression of individual genes in blood of young healthy pigs were identified for concurrent growth rate and subjective health scores prior to the challenge, and for mortality, a combined mortality-treatment trait, and feed conversion rate after the challenge. Gene set enrichment analyses revealed three groups of gene ontology biological process terms that were related to disease resilience: 1) immune and stress response-related terms were enriched among genes whose increased expression was unfavorably associated with both pre- and post-challenge traits, 2) heme-related terms were enriched among genes that had favorable associations with both pre- and post-challenge traits, and 3) terms related to protein localization and viral gene expression were enriched among genes that were associated with reduced performance and health traits after but not before the challenge. CONCLUSIONS: Gene expression profiles in blood from young healthy piglets provide insight into their performance when exposed to disease and other stressors. The expression of genes involved in stress response, heme metabolism, and baseline expression of host genes related to virus propagation were found to be associated with host response to disease.

Reference Transcriptomes of Porcine Peripheral Immune Cells Created Through Bulk and Single-Cell RNA Sequencing.

Pigs are a valuable human biomedical model and an important protein source supporting global food security. The transcriptomes of peripheral blood immune cells in pigs were defined at the bulk cell-type and single cell levels. First, eight cell types were isolated in bulk from peripheral blood mononuclear cells (PBMCs) by cell sorting, representing Myeloid, NK cells and specific populations of T and B-cells. Transcriptomes for each bulk population of cells were generated by RNA-seq with 10,974 expressed genes detected. Pairwise comparisons between cell types revealed specific expression, while enrichment analysis identified 1,885 to 3,591 significantly enriched genes across all 8 cell types. Gene Ontology analysis for the top 25% of significantly enriched genes (SEG) showed high enrichment of biological processes related to the nature of each cell type. Comparison of gene expression indicated highly significant correlations between pig cells and corresponding human PBMC bulk RNA-seq data available in Haemopedia. Second, higher resolution of distinct cell populations was obtained by single-cell RNA-sequencing (scRNA-seq) of PBMC. Seven PBMC samples were partitioned and sequenced that produced 28,810 single cell transcriptomes distributed across 36 clusters and classified into 13 general cell types including plasmacytoid dendritic cells (DC), conventional DCs, monocytes, B-cell, conventional CD4 and CD8 αß T-cells, NK cells, and γδ T-cells. Signature gene sets from the human Haemopedia data were assessed for relative enrichment in genes expressed in pig cells and integration of pig scRNA-seq with a public human scRNA-seq dataset provided further validation for similarity between human and pig data. The sorted porcine bulk RNAseq dataset informed classification of scRNA-seq PBMC populations; specifically, an integration of the datasets showed that the pig bulk RNAseq data helped define the CD4CD8 double-positive T-cell populations in the scRNA-seq data. Overall, the data provides deep and well-validated transcriptomic data from sorted PBMC populations and the first single-cell transcriptomic data for porcine PBMCs. This resource will be invaluable for annotation of pig genes controlling immunogenetic traits as part of the porcine Functional Annotation of Animal Genomes (FAANG) project, as well as further study of, and development of new reagents for, porcine immunology.

Effect of ARTEMIS (DCLRE1C) deficiency and microinjection timing on editing efficiency during somatic cell nuclear transfer and in vitro fertilization using the CRISPR/Cas9 system.

The ability to efficiently introduce site-specific genetic modifications to the mammalian genome has been dramatically improved with the use of the CRISPR/Cas9 system. CRISPR/Cas9 is a powerful tool used to generate genetic modifications by causing double-strand breaks (DSBs) in DNA. Artemis (ART; also known as DCLRE1C), is a nuclear protein and is essential for DSB end joining in DNA repair via the canonical non-homologous end joining (c-NHEJ) pathway. In this work, we tested whether ART deficiency affects DNA repair following CRISPR/Cas9 induced DSBs in somatic cells. We also demonstrated the effect of microinjection timing on embryo developmental ability and gene targeting efficiency of CRISPR/Cas9 system to disrupt the interleukin 2 receptor subunit gamma (IL2RG) locus using porcine in vitro fertilization (IVF) and somatic cell nuclear transfer (SCNT) derived embryos. In comparison to non-injected controls, CRISPR/Cas9 injection of IVF derived zygotes at 4 h and 8 h after fertilization did not impact cleavage and blastocyst rate. Gene modification rate was observed to be higher, 53.3% (9/16) in blastocysts injected 4 h post-fertilization as compared to 11.1% (1/9) in blastocysts injected 8 h post-fertilization. Microinjection 8 h after chemical activation of SCNT derived embryos decreased blastocyst development rate compared to non-injected controls but showed a higher gene modification efficiency of 66.7% as compared to 25% in the 4 h post-activation injection group. Furthermore, we observed that male ART-/- and ART+/- porcine fetal fibroblast (pFF) cells showed lower modification rates (2.5% and 1.9%, respectively) as compared to the ART intact cell line (8.3%). Interestingly, the female ART-/- and ART+/- pFF cells had modification rates (4.2% and 10.1%, respectively) similar to those seen in the ART intact cells. This study demonstrates the complex effect of various parameters such as microinjection timing and ART deficiency on gene editing efficiency in in vitro derived porcine embryos.

Gene expression in tonsils in swine following infection with porcine reproductive and respiratory syndrome virus.

BACKGROUND: Porcine reproductive and respiratory syndrome (PRRS) is a threat to pig production worldwide. Our objective was to understand mechanisms of persistence of PRRS virus (PRRSV) in tonsil. Transcriptome data from tonsil samples collected at 42 days post infection (dpi) were generated by RNA-seq and NanoString on 51 pigs that were selected to contrast the two PRRSV isolates used, NVSL and KS06, high and low tonsil viral level at 42 dpi, and the favorable and unfavorable genotypes at a genetic marker (WUR) for the putative PRRSV resistance gene GBP5. RESULTS: The number of differentially expressed genes (DEGs) differed markedly between models with and without accounting for cell-type enrichments (CE) in the samples that were predicted from the RNA-seq data. This indicates that differences in cell composition in tissues that consist of multiple cell types, such as tonsil, can have a large impact on observed differences in gene expression. Based on both the NanoString and the RNA-seq data, KS06-infected pigs showed greater activation, or less inhibition, of immune response in tonsils at 42 dpi than NVSL-infected pigs, with and without accounting for CE. This suggests that the NVSL virus may be better than the KS06 virus at evading host immune response and persists in tonsils by weakening, or preventing, host immune responses. Pigs with high viral levels showed larger CE of immune cells than low viral level pigs, potentially to trigger stronger immune responses. Presence of high tonsil virus was associated with a stronger immune response, especially innate immune response through interferon signaling, but these differences were not significant when accounting for CE. Genotype at WUR was associated with different effects on immune response in tonsils of pigs during the persistence stage, depending on viral isolate and tonsil viral level. CONCLUSIONS: Results of this study provide insights into the effects of PRRSV isolate, tonsil viral level, and WUR genotype on host immune response and into potential mechanisms of PRRSV persistence in tonsils that could be targeted to improve strategies to reduce viral rebreaks. Finally, to understand transcriptome responses in tissues that consist of multiple cell types, it is important to consider differences in cell composition.

Differential induction of innate memory in porcine monocytes by ß-glucan or bacillus Calmette-Guerin.

Innate immunomodulation via induction of innate memory is one mechanism to alter the host's innate immune response to reduce or prevent disease. Microbial products modulate innate responses with immediate and lasting effects. Innate memory is characterized by enhanced (training) or depressed (tolerance) innate immune responses, including pro-inflammatory cytokine production, to secondary exposure following a priming event. To investigate the ability of ß-glucans and bacillus Calmette-Guerin to induce innate training or tolerance in pig cells, porcine monocytes were cultured with priming agonist (ß-glucans or bacillus Calmette-Guerin) then re-stimulated 5 d later with a heterologous microbial agonist to determine induction of innate memory. Priming with ß-glucan from Saccharomyces cerevisiae depressed IL-1ß and TNF-α cytokine responses to re-stimulation with LPS, indicative of a tolerized state. However, bacillus Calmette-Guerin priming induced a trained state in porcine monocytes, as LPS re-stimulation enhanced IL-1ß and TNF-α gene expression and protein production. We present the first evidence of innate memory in pig monocytes, with bacillus Calmette-Guerin (training) or Saccharomyces cerevisiae ß-glucan (tolerance). Induction of a trained or tolerized state in vitro is a first step to identify agonists to alter the innate immune system at the animal level with the intent of enhancing disease resistance.

A comparative analysis of chromatin accessibility in cattle, pig, and mouse tissues.

BACKGROUND: Although considerable progress has been made towards annotating the noncoding portion of the human and mouse genomes, regulatory elements in other species, such as livestock, remain poorly characterized. This lack of functional annotation poses a substantial roadblock to agricultural research and diminishes the value of these species as model organisms. As active regulatory elements are typically characterized by chromatin accessibility, we implemented the Assay for Transposase Accessible Chromatin (ATAC-seq) to annotate and characterize regulatory elements in pigs and cattle, given a set of eight adult tissues. RESULTS: Overall, 306,304 and 273,594 active regulatory elements were identified in pig and cattle, respectively. 71,478 porcine and 47,454 bovine regulatory elements were highly tissue-specific and were correspondingly enriched for binding motifs of known tissue-specific transcription factors. However, in every tissue the most prevalent accessible motif corresponded to the insulator CTCF, suggesting pervasive involvement in 3-D chromatin organization. Taking advantage of a similar dataset in mouse, open chromatin in pig, cattle, and mice were compared, revealing that the conservation of regulatory elements, in terms of sequence identity and accessibility, was consistent with evolutionary distance; whereas pig and cattle shared about 20% of accessible sites, mice and ungulates only had about 10% of accessible sites in common. Furthermore, conservation of accessibility was more prevalent at promoters than at intergenic regions. CONCLUSIONS: The lack of conserved accessibility at distal elements is consistent with rapid evolution of enhancers, and further emphasizes the need to annotate regulatory elements in individual species, rather than inferring elements based on homology. This atlas of chromatin accessibility in cattle and pig constitutes a substantial step towards annotating livestock genomes and dissecting the regulatory link between genome and phenome.

Changes in H3K27ac at Gene Regulatory Regions in Porcine Alveolar Macrophages Following LPS or PolyIC Exposure.

Changes in chromatin structure, especially in histone modifications (HMs), linked with chromatin accessibility for transcription machinery, are considered to play significant roles in transcriptional regulation. Alveolar macrophages (AM) are important immune cells for protection against pulmonary pathogens, and must readily respond to bacteria and viruses that enter the airways. Mechanism(s) controlling AM innate response to different pathogen-associated molecular patterns (PAMPs) are not well defined in pigs. By combining RNA sequencing (RNA-seq) with chromatin immunoprecipitation and sequencing (ChIP-seq) for four histone marks (H3K4me3, H3K4me1, H3K27ac and H3K27me3), we established a chromatin state map for AM stimulated with two different PAMPs, lipopolysaccharide (LPS) and Poly(I:C), and investigated the potential effect of identified histone modifications on transcription factor binding motif (TFBM) prediction and RNA abundance changes in these AM. The integrative analysis suggests that the differential gene expression between non-stimulated and stimulated AM is significantly associated with changes in the H3K27ac level at active regulatory regions. Although global changes in chromatin states were minor after stimulation, we detected chromatin state changes for differentially expressed genes involved in the TLR4, TLR3 and RIG-I signaling pathways. We found that regions marked by H3K27ac genome-wide were enriched for TFBMs of TF that are involved in the inflammatory response. We further documented that TF whose expression was induced by these stimuli had TFBMs enriched within H3K27ac-marked regions whose chromatin state changed by these same stimuli. Given that the dramatic transcriptomic changes and minor chromatin state changes occurred in response to both stimuli, we conclude that regulatory elements (i.e. active promoters) that contain transcription factor binding motifs were already active/poised in AM for immediate inflammatory response to PAMPs. In summary, our data provides the first chromatin state map of porcine AM in response to bacterial and viral PAMPs, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the role of HMs, especially H3K27ac, in regulating transcription in AM in response to LPS and Poly(I:C).